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O ' Abstract 

^ : 

■ A rigorous treatment is given of the Green's function of the = 1 supersym- 

metric heat equation in one spatial dimension with a distribution initial value. The 
asymptotic expansion of the supersymmetric Green's function as t tends to 0^ is 
also derived. The coefficients of the expansion generate all the members of the 
r~| , supersymmetric = 1 KdV hierarchy. 

> : 

^ : 1 Introduction 



We analize in this work the quantum mechanical problem, in the "euclidean" domain, 
associated with supersymmetric integrable models. The N = 1,2 supersymmetric exten- 
sions of the KdV equation were found several years ago Pfl |21 El El- 
Later, the = 3,4 cases were also discussed ISHEIIZI- 

A lot of work has been done on the subject, in particular the bihamiltonian struc- 
ture of the super KdV equations was studied in jHl IHl- For a review of the = 1,2 
supersymmetric KdV equations see [T^. Extensions of the supersymmetric model were 
proposed in ^T]. Supersymmetric quantum problems in other contexts are also relevant. 
In particular, the quantum mechanical problem arising from the regularization of the 
Supermembrane in eleven dimensions. The nature of the spectrum of the corresponding 
hamiltonian has been described in fT^ IT^ ITU [T3] , however a satisfactory understanding 
of the corresponding Green's function has not been achieved. 
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In ^ni; using the exact sequence for the = 1 supersymmetric ring, exphcit expres- 
sions for the Green's function of the = 1 supersymmetric heat equation were obtained 
and using a formal asymptotic expansion of it, when t ^ 0^, all the supersymmetric 
N = 1 KdV hierarchy was generated. In the present work we prove the convergence of 
the series defining the supersymmetric Green's function with a distribution initial value 
and give a rigorous treatment of its asymptotic expansion as t — 0"*". 

We thus complement our previous work with rigorous statements. In section 2 we 
present general considerations, the new results are developed in sections 3 and 4. 

2 The Susy Heat Operator 

The usual homogeneous heat equation with potential may be written 

Hu = (j)u 

with H = dt — and u{t,x), real-valued functions of real variables. In the su- 
persymmetric extension to be considered here, the functions will take their values in an 
exterior algebra A. In more detail: 

Let A be generated by anticommuting elements ei, . . . , e.^ which satisfy = 0. The 
subspaces of bosonic and fermionic elements of A are the ±1 eigenspaces of the involution 
A — >• A wich sends the generators Cj to — Cj. 

For 1 < i < m the operator (9j : A — > A is defined by the equations diX = and 
A = CjCt -|- P, the elements a and P being independent of the generator Cj. This operator is 
a super derivation because it anticommutes with the involution and satisfies the product 
rule di{XiJ,) = {diX)fi + X{difi). 

Let Ai C A be the subalgebra generated by ei,...e„ where n < m, and let 9 = 
{e„+i, . . . , Cm} be the set of remaining generators. If 6' G 6 then the superderivation 
dg : A A also acts in the ring of all functions / : M ^ A. 

If the elements of G are written as 61,62, .. . then the functions f : MP ^ A can be 
written in the notation 

f{x,6) = J2(^,,---6,Jj{x) 

with J = (ji < ■ ■ ■ < jg) and fj:W^Ai. One sees that dgf has the same fj but with 
different coefficients. Evidently 6 could have any size; but in what follows it will have no 
more than three elements 6, 6', 6". 

Supposing now that p = l,n + 1 = m and = {6} we begin by observing that the 
ring of functions / : R ^ A admits the operator 

D = de + 6d,, 

at least if we restrict to the case f{x,6) = fi{x) + 6f2{x) with different iable /i and /2. 
Then Df = f2{x) + 6^J^{x). 
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The functions u{t,x) and (f){x) in the usual heat equation are now replaced by 

u{t, X, 0) — Ui{t, x) + 9u2{t, x) 

(t){x, e) = a{x) + e(3{x) 

with Ml, M2, a, and /? taking their values in Ai. The SUSY heat equation, with the same 
H — dt — dl as before, is defined to be 

Hu = D{(j)u). 

In terms of the components it is 

Hui = Pui + au2 
Hu2 — dx{a.ui). 

The usual heat equation Hui = j3ui reappears if a = and the values of j3 and ui 
are multiples of the identity clement of A. This requires j3 to be of even type, that is, 
j3 = 13. This possibility is admitted by assuming that the SUSY potential function is 
fermionic, that is, = —0. Then a = —a, and the componentwise version of the SUSY 
heat equation is 

Hui — (5ui — au2 
Hu2 = d:c{dui). 

Going back for a moment to the usual heat equation we recall that the inhomogenous 
equation 

Hu — (f)u-\- f, 

with /(t, x) defined in the upper half plane i > 0, is solved by 

l*t /"OO 

u{t,x)^ / g{t - s,x,y)f{s,y)dyds 

Jo J -oo 

where g{t, x, y) is the Green's function of the equation: it satisfies Hg — (pg as a. function 
of X and t when t > 0, and has the initial value 

\im g{t,x,y) = 5(x-y), 

the distributional Dirac delta function. 

An analogous SUSY Green's function G is now to be developed for solving the SUSY 
equation Hu — D{(j)u) + Df. There are three steps: first a SUSY integration formula is 
developed for solving Hu — Df; second an exphcit construction is given for the SUSY 
Green's function G with attention to its symmetries; and third it is shown that G has an 
asymptotic expansion in powers of t. 
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Guided by the usual g{t, x, y) we might suggest that the SUSY Green's function should 
depend (x, 6) and (y, 9') in place of x and y. Accordingly 

G{t, X, y; e, e') = a + eb + O'c + 99' d 

where a, 6, c, and d are Ai-valued functions of x and y, Then, given f{t,x,9), an 
analogous definition of a;, 9) would be 

u{t,x,9)= / G{t- s,x,y;9,9')f{s,y,9')d9'dyds. 

Jo J-oo 

The integration with respect to 9' is to be understood as the partial derivative with respect 
to 9': that is, given in general any J{xi, . . . ,Xp;9i, . . . , 9n), its integral with respect to 9j 
is defined to be 

With f{t, X, 9) — fi{t, x) + 9f2{t, x), the integration d9' leaves the integrand in the form 

{ch + af2)+9{-dh-bf2). 

Therefore, with u{t, x, 9) — Ui{t, x) + 9u2{t, x) we have 

ui = Jc/i + / a/2 
U2^-Jdfi- J 6/2 

with J signifying the remaining integrations dyds. 

Since Df = /2 + 9dxfi, the equation Hu = Df will be satisfied if we choose b — c — 0, 
while a and d should satisfy 

H{Jaf2) = f2 
H{Jdh)^-dJ,. 

Therefore a should be the standard Gauss bell heat kernel and d should be the negative 
of its derivative with respect to x. With the definition 

Go{t, x-y;9, 9') = ^e"^ (l + 99'^) 

we have a SUSY Green's function which solves Hu = Df hj u = J G^f . 

The defining properties of Go are that it should satisfy HGq = for t > 0, with initial 
value 

limGo = -y)- 00' 5' {x - y), 
with S' being the distributional derivative of the Dirac delta function. 
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In addition one observes that Go is bosonic (Go = Go) and enjoys the symmetry 

Go{t,x,y;e,e')=Go{t,y,x;e',e). 
Turning now to the equation 

Hu = D{(t)u) + Df 

with a fermionic potential <^{x,0) — a{x) + 6(3{x), we ask for its Green's function. Fol- 
lowing what was seen for Gq, we should ask that G{t,x,y;6,6') satisfy HG = D{(j)G) 
as a function of t, x, and 6 when t > 0, and that it should have initial value limG = 

m 

S{x — y) — 99' 6' {x — y). Supposing that such a G exists, we argue that u = J Gf satisfies 
Hu = D{(t)u) + Df. 
Writing 

G^ A + OB + e'c + ee'D 

with A, B, C, and D functions of t, x, and y, we note that the componentwise version of 
HG = D{(I)G) is 

HA = pA- aB 
HB = d^^{aA) 
HC = aD + I3C 
HD^-d^{aC). 

(Rewritten with explicit independent variables the first equation, for example, is 

{dt - dl)A{t, X, y) = (5{x)A{t, X, y) - a{x)B{t, x, y).) 
The parities A = A, D = D, S = -S and G = -G give 

«l = /G/i + /A/2 

U2 = jBh-jDh. 
To apply H we use the general fact that 

"t POD 



v{t,x)^ / K{t - s,x,y)e{s,y)dyds 

Jo J-oo 

implies 



pt POQ POO 

Hv{t,x)= / {HK){t-s,x,y)e{s,y)dyds + lim K{s,x,y)e{t,y)dy. 

Jo J-oo J-oo 



'0 J-oo 

The initial values imposed on G leave 



Hm^f2 + JiHC)h + f{HA)f2 
Hu2 = dj, + J{HB)f2 - J{HD)f,. 



However 

J{HA)f2 = Jo JZ, iBix)A{t - s, X, y)h{s, y) - a{x)B{t - s, x, y)h{s, y))dyds 

= (3jAf2-ajBf2. 

In like manner 

J{HB)f2 = /o X!^ d^a{x)A{t - s, x, y)/2(s, y)dyds 
= d^aj Af2. 

Similar computations for HC and HD lead to 

Hui = /2 + (3ui - au2 
Hu2 = d^fi + d:,{aui), 

which is just the componentwise version of Hu = D{4>u) + Df, as desired. 



3 Construction of the SUSY Green's Function 

A potential (f){x, 9) — a{x) + 6(5{x) is given. In addition to the conditions a — —a and 
f3 = [3, we assume that a and f3 are bounded infinitely differentiable functions from M to 
A, and that all their derivatives are bounded. 
Formally we put 

oo 

G = Gm 

m=0 

where Gq, defined in the preceding section, satisfies HGq = in the upper half plane and 
has the initial value 5{x — y) — 99'S'{x — y). 

The desired equation HG — D{(f)G) suggests the recursion -ffm+i = D{(f)Gm)- This 
equation in turn is solved by integration with Go, as follows. 

For fixed y, 9' the function (pG^ takes the form 

f{t,x,9)^<l>ix,9)Gm{t,x,y;9,9') 

for variable x, t, 9. To solve Hu = Df we replace x, t and 9 by variables of integration 
u, s and 9". Then 

u{t, x,9)^ / Go{t -s,x-u;9, 9")f{s, u, 9")d9"duds, 

Jo J-oo 

giving 

ft fOO 

Gm+i{t,x,u;9,9') = / / Go{t — s,x — u;9,9")(l){u,9")Gm{s,u,y;9",9')d9"duds. 

Jo J-oo 
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We may write 

Gm — Ajji + OBjn + O'C^ + 99'Djn 
in general, where for m = we already know that Gq — g — 99' h, with g{t,x — y) — 

(x—y)'^ o 

72^e ^ and h{t,x-y) = £g{t,x-y). 

The formal integration d9" is done by applying ^ to 

{g - 99"h){a + 9"P){A^ + 9" + 9' Cm + 9" 9' Dm). 

This leaves 

{gpAm - gaBm) + 9{haAm) + 9\gl3Cm + - 99\haCm). 

Therefore except for the remaining integrations duds we have the componentwise recursion 
formula: 

Am+l = gf^Am - gOtBm 

Bm+\ — ha Am 

Cm+1 = gPCm + gaDm 

Dm+1 — —haCm- 
This can be recast in matrix form as 

Am+l -Cm+l \ ^ f 9/3 -gOi \f Am -Cm 
Bm+1 Dm+1 J \ha J \ Bm Dm 

giving at least formally the explicit formula 

Am -Cm \ ^ ( g(3 -g(^ VV ^ ^ 

Bm Dm J \ha J {O -h 

The remaining integrations of the real variables u and s must now be examined. The 
equation Bi = hag for example means that 



Bi{t,x,y)= I I h{t — s,x — u)a{u)g{s,u — y)duds. 

Jo J-oo 

The general case involves various products of the four functions a, l3, g, and h. To 
treat it one can define a "word" to be any clement of a finite Cartesian product 
ExflxExflx---xflxE where E — {g, h} and fl — {a, /?}. Then, given 

Wm = eoUJieiU2 ■ ■ ■UJme-m, 

a function oit, x and y arises through the integral 



Wm{t,X,y) = / eo(t - Si,X - Ui)ui{ui)ei{si - S2U1 - U2) ■ ■ ■ UJm{Um)em{Sm,Um " V) 
'Im{t) 



where Im{t) C M^"* is the set of (si, ui, S2, U2, ■ ■ ■ , Sm, Um) satisfying t > si > S2 > • • • > 
Sjji > 0. The convergence of this integral and of the series defining G are now to be 
discussed, as well as the distribution initial value of G as t tends to zero. 

First we observe that g and h are in the space JC of real-valued continuous functions 
k{x, t) in the upper half plane which satisfy, for all < Ti < T2, 

lo' I-oo\Hs,u)\duds < 00 
sup \k{s, u)\ < 00. 
—00 < u < 00 
Ti < s < Ta 

The space /C is a commutative ring under convolution 

{ki * k2){t,x) — / / ki{t — s,x — u)k2{s,u)duds. 

Jo J-00 

If we take k{t, x) to be some positive multiple of g{t, x) + \ h{t, x)\, we obtain the estimate 

\Wm{t,x,y)\ < l{t,x-y) 
in which / is the (m + l)-fold convolution 

This shows that Wm{t,x,y) is in /C for any value of y, and is defined by an absolutely 
convergent integral. 

In order to estimate I in points of the strip < i < 1 we define the three auxiliary 
functions 

1*00 

ki{t) = / / k{s,u)duds 

Jo J-00 

T{t) = sup k{s, u) 
—00 < u < 00 

t < S<1 

p(r) = sup k{s, u) 
u < —r or u > r 
< s < 1 

Then from the definition of the (m + l)-fold convolution one can deduce 



l{x,t) < (m + l)(A;i(t))"V 
l{x,t)<{m+l){k,{t)rp(^). 
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For some C the function k{t,x) satisfies 

C 

k(t,x) < — 

2 

in the strip < t < 1, as can be seen by using e'' > 1 + r + y in the formulas for g and h. 
Hence r(t) < f and p(r) < Also ki{t) < Cy/i, again by checking the same formulas. 
This gives us the estimates 

l{x,t) < {m + lfC'^+H'^ 

m 

l{x,t) < {m + lfC^+^i^. 

oo 

The first estimate gives the absolute convergence of any series Wmjt, x, y) in some 

1 

oo 

strip, say < i < Taking off the first two terms we also see that limy~] Wm — 0, 
uniformly in x and y. 



m=3 



oo 

In order to show that the full series satisfies hm > Wm = in the distribution sense 

m=l 

we must examine Wi and W2- 

The second estimate shows, for any m > 1, that limWm{t,x,y) — uniformly in any 

oo 

set of X and y of the form \x — y\ > r > 0. Therefore the desired initial value of Wm 

m=l 

would follow if Wi and W2 were bounded functions. 

But this is clear for W2, by the first estimate. For Wi wc note that the recursion 
formula for Gi only mentions g^g, gah, and hag; therefore \h\ * \h\ need not be analyzed. 

One can easily check 



{lZo\9is^u)\^duj ^Cis 5 

(rjh{s,u)rduy" = c2s-i 

for some constants, giving 



J — c 



g{t — s,x — u)\h{s,u\du < CiC2{t — s) 45 4 



by the Schwartz inequality. Then because J^{t — s)~^s~^ds is a finite positive constant 
independent of t > 0, we conclude that g * \h\ is a bounded function for —00 < x < 00, 
t > 0. (The remaining case g * g — tg{x, t) is obvious.) 
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To sum up: / ^ Wmjt, x, y) is absolutely and uniformly convergent in some strip < 
1 

t < e, and has distribution limit zero as t tends to zero. 

Prom the recursion formulas for the Gm it is clear that each of A^, • • • , -Dm is a linear 
combination of no more than 2™ words Wm, with ±1 coefficients. Therefore the SUSY 
Green's function is given by the series 

oo 

G{t,x,y;e,e') = J2Gm, 



convergent in a smaller strip < t < e/2, with the same initial value as the first term Gq: 

limG = 6{x-y)-ee'6'{x-y) 

The construction of G being complete, we turn now to its symmetries. These depend 
on the words produced by the recursion algorithm. If Wm — eQUJ\e\ijJi ■ ■ •cum^rn- Then 

Wm{t,x,y)^±W:^{t,y,x) 

where 

= emUJm • ■■uj2eiUJieo 

is the word Wm written backwards, and the sign depends on how many times h and a 
appear in the word. Here is the proof: 

For fixed t > the involution of R given by s — > s = t—s produces a volume-preserving 
self-map of the set Irn {t) C IR^™ where the integration takes place: (si, Ui, S2, U2, ■ ■ ■ , SmUm) 
is sent to (s^, Mm; • • • , Si, Ui). After this change of variable the new integrand is the prod- 
uct Oi eQ{t -Sm, X - Um)ei{Sm-Sm-l,Um- Um-l) ■ ' ' emisi, Ui - y) and LL!i{Um) ■ ■■LVmiUi). 

At this point we note that g and h satisfy g{si — S2;ui — U2) — g{s2 — si,U2 — ui) 
and h{si — S2,Ui — U2) = —h{s2 — 'si,U2 — Ui) when t > Si > S2 > 0. Consequently the 
products of g^s and /I's before and after the change of variable differ by the sign (—1)'^, n 
being the number of times h appears among {cq, Ci, . . . , 6^}- 

On the other hand the products of a's and /3's satisfy 

i^(i^-i) 

UJi{Un) ■ ■ ■ UJn{ui) ^ {-1) 2 LUn{ui) ■ ■ ■ UJi{Un) , 

V being the number of times the fermionic term a appears among {a;i, . . . , a;^}- The 
factors Cfc commute with all elements of the exterior algebra, so after one more reordering 
we obtain 

Wm{t,x,y) = W*Jt,y,x)i-in-iy-^. 

The quantity e{W) = (— 1)^(— 1) * 2 ' is well-defined for any word in the four-letter 
alphabet {g,h,a,P} and satisfies 

e{WiW2) = e{Wi)e{W2){-iy"''. 
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The recursion relation 



gives 



and 



Bm+1 Dm+1 J V ha Am -haCm 



Ai -Ci \ ^ f gl3g gah 
B, Di J I hag 



A2 -C2 \ ^ f gPgPg - gahag gpgah 
B2 D2 J \ hagPg hagah 



Inspection shows that e{Am) ~ +1 while e{Bm) = ^(C'm) = ^(-0^) = — 1 for m = 1 
and 2, and that A^ and are bosonic while B^, and are fermionic, this referring to 
the parity of the number of times v that the letter a appears in the word. 

The recursion relation can be used to show that these statements hold for all m, 
namely, that B^, and Cm are formal linear combinations of fermionic words satisfying 
£ = — 1, while Am and Dm are combinations of bosonic words with £ = 1, resp. e — —1, 
for the terms in the two cases. (For example e{ha) — —1, e{gah) — —1 giving e{hagah) — 
—e{ha)s{gah) — —1 because both factors are fermionic.) 

The examples suggest that = Am, = Dm, and B^ = —Cm- To see that this is 
true for all n one can write the recursion formula in expanded form as 

Am -Cm \ ^( g \ f P -a \ f g 0\ f P -a\(g 
Bm Dm J \0 -h J \-a J \ -h J "' \ -a J \ -h 

which gives 

Am+2 -Cm+2 \ ^ ( dP "fi'^ \ f Am 'Cm \ f jSg ah 
, Bm+2 Dm+2 J \ ha J \ Bm Dm J \ -ag 

and hence 

Am+2 = gf3Am(3g - gaBmPg + gf^CmOtg + gaDmag 
Bm+2 = haAm(5g + haCmOtg 
Cm+2 = gaBmOth - g(3Amah 
Dm+2 = haAmah. 

The desired conclusion A^_^2 — ■Am+2, — Dm+2, — —Cm+2 now follows from 

the same statement for m. 

Combining the above we obtain 

Am{t,x,y) = Am{t,y,x) 
Bm(t,x,y) = Cm(t,y,x) 
Dm{t,x,y) = -Dm{t,y,x). 
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With Gm{t, X, y, 9, 9') = Am + 9Bm + 9' Cm + 99' Dm these three equations take the form 

Gm{t,x,y,9,9')^Gm{t,y,x,9',9). 

oo 

The same condition for the SUSY Green's function G = ^ Gm is the expression of its 



symmetry with respect to x y and 9 9' interchanges. 

4 The Asymptotic SUSY Green's Function 

Computations of asymptotic series will be done in the following setting. 

Let fl be an open subset of M". Then the space C^{Q) is to consist of the infinitely 
differentiable functions defined in fl which, together with all partial derivatives, admit 
continuous extensions to the closure of Q and are bounded by polynomials. 

When f2 C is the half-space t > 0, L E C^{fl) admits the asymptotic expansion 

oo 

L{t,x,y) ~ ^rL„(x,t/) 

n=0 

in which Ln{x, y) — ^ \im d^L{t, x, y). The error estimate is 

N 

L{t, x,y) = J2 t"Ln{x, y) + x, y) 

n=0 

where, by Taylor's Theorem Rn+i G Cj^(f2). 

Given K G C|^(]R^), the square root substitution may produce an element of C'^{VL): 
if K is defined everywhere in and satisfies K{t,x,y) = K{—t,x,y) then K{\/t, x,y) e 

c^m. 

When elements of these spaces appear in improper integrals the convergence will be 
deduced from the polynomial growth conditions. These permit separation of variables, in 
that every polynomial satisfies 

\P{x„...,x^)\<C(l+xf)---(l + xf) 

for some C > 0, A: > 0. 

These preparations done, we examine the transformation L ^ M given by 



g{t,x -y)M{t,x,y) = 




g{t — s,x — u)L{s, u, y)g{s, u — y)duds. 
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The Gauss bell heat kernel admits the factorization 

g{t - s,x- u)g{s, u-y) = g{t, x - y)g{T, u - w) 

in which < s < t, s — at, t — a{l — a)t, and w — ax + {1 — a)y. Thus, except for a 
constant factor, 

M{t, X, y) = L{s, u, y)g{T, u - w)duds 

= Jo IZo ^('^' ^ + v^^' y)e'^dvds 

after the change of variable u = w + ^/rv, du = ^/Tdv. 

By hypothesis L{t,x,y) G (H). Then L(s, w + Tv,y) e C^(Ai), where A C is 
defined by s > 0. The polynomial bound gives the absolute convergence of 

J{s,T,w,y) ^ / L{s,'w + Tv,y)e ^ dv, 

J —oo 

a function in C^{A2), with A2 C defined by s > 0. 

Evidently J is invariant under r — r. Therefore J(s, ^,w,y) = K{s,T,w,y) is in 
C|°(A3), with A3 C given by s > 0, r > 0. This leaves M in the form 

M{t, X, y) = /p K{s, T, w, y)ds 
= i K{at, a{l - a)t, ax + {1 - a)y, y)da, 

after the change of variable s ~ at, ds — tda. 

It is now evident that the construction L ^ M is a linear transformation into itself 
of the vector space Cj^(fi), with C M'^ given by t > as before. Furthermore if 

L=Y, t'^Lnix, y) then M = rM„(x, y). 

n>N n>N+l 

When dx or dy is applied to the equation at the very beginning, the integrand gLg is 
replaced by hLg or gLh, and M is replaced by dxM ± ^^M, for example. This function 
is in C^IQ) once more, but the cancellation of t means that the resulting asymptotic 

series may, like that of L, involve all with n > N. 

00 

The Green's function G — Gm can now be expanded term by term. In general, 

m=0 

words W — eQUJieiUJ'2, ■ ■ -00^^^ give functions of three variables W{t, x, y) by the multiple 
integral formula. If W is lengthened to eujW with e — g or h and uj — a or (5, one has 
the recursion formula 



{eujW){t,x,y) = 




e{t — s,x — u)u{u)W{s, u, y)duds. 



If W{t,x,y) — g{t,x — y)L{t,x,y) for some L e C^(Q) then {eu}W){t,x,y) — g{t,x — 
y)M{t, x.y) where M G C^{Q) is the result of applying the L ^ M transformation to 
u!{x)L{t,x,y) and then taking a derivative if necessary. 



13 



The observations about lowest powers of t permit the specific formula 

oo 

W{t,x,y) = g{t,x-y) ^""^nix^y) 

n=N-l 

where = N{W) is the number of times that g occurs in the word W. Therefore, by 
induction, every word gives rise to an asymptotic expansion as written, with coefficients 
L„(x,y)eCf(R2). 

Then, looking at Gi, we see that 

oo 

G^(t, X, y, e, e') = g{t, x - y)Y, f^L^^ix, y, 6, 9'). 

n=0 

Each Ln is a linear combination of four A-valued elements of C|^(M^), with coefficients 

1, 9, 9', and 99'. In particular Lo(a;, |/, 0, 0) = 0, because, with Ai — gPg, we have 
Ai(0,x,y) = 0. 

When we substitute 

g{t, x-y)^ Go{t, x-y,9, 9') (l - ^^'^\~^^ ) , 

no negative powers of t appear, because 99'Lo{x, y, 9, 9') — 99'Ai{Q., x, y) = 0. Therefore 

oo 

G,{t, X, y, 9, 9') = G'o(i, x-y,9, 9') ^ rL„(x, y, 9, 9'), 

n=0 

with coefficients in C|^(M^) as before. 

In the general case of Gm, the four entries in the matrix recursion formula are finite 
linear combinations of words with integer coefficients. For such an entry, define N to be 
the minimum of N{W) for the words appearing in the entry. 

Suppose for some m that N{Am) > p + 1 and N{Bm), N{Cm), N{Dm) > P for some 
nonnegative integer p. This holds for m = 1 and m — 2 with p — 1, by inspection of 
the formulas. An examination of the m ^ m + 2 recursion formula then shows that 
N{A^+2) >P + 2 while Ar(S„+2), N{G^+2), N{Dm+2) >P+l. 

This shows that 

lim N{Gm) = +00, 

m— »oo 

with N{Gm) being the minimum of N{A^), . . . , N{D^). Thus, with a change of notation, 
Gm{t, X, y, 9, 9') = Go{t, x-y,9, 9') ^ rL^,„(x, t/, 9, 9'). 
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Passing to we note that only finitely many m can make a contribution to any 

m=0 

given power of t. Therefore 

oo 

G(t, X, y, e, 6') ~ Go(t, x-y,e, 6') f'M^ix, y, 6, 6'), 

n=0 

each Mn being a combination of four elements of C|°(]R^). For example 

y QQ' ry 

a{s)ds + — / P{s)ds. 

The coefficients Mn can be determined from a recursive system of differential equations 
derived from the G = t^Mn equation and the differential equation 

EG = D{(j)G). 

This completes the verification that the SUSY Green's function possesses an asymptotic 
expansion, in the above sense. 

In it was shown that this asymptotic expansion, which is now rigorously estab- 
lished, generated all the members of the = 1 supersymmetric KdV hierarchy. 

5 Conclusions 

A rigorous treatment has been given of the Green's function G of the = 1 super- 
symmetric heat equation in one spatial dimension with a distribution initial value. The 
convergence of the series defining G is proven as well as the existence of its asymptotic 
expansion as t 0^. It was shown in [16J that the formal asymptotic expansion generates 
all the members of the supersymmetric = 1 KdV hierarchy. The asymptotic expansion 
for the Green's function of the KdV equation was discussed in |17j . 
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